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Role of the coin in the spectrum of quantum walks 
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Abstract 

The most elementary quantum walk is characterized by a 2-dimensional unitary coin 
flip matrix, which can be parameterized by 4 real variables. The influence of the choice 
of the coin flip matrix on the time evolution operator is analysed in a systematic way. 
By changing the coin parameters, the dispersion and asymmetry of eigenvalues of the 
time evolution operator can be tuned in a controlled way. The reduced eigenvectors in 
coin space are distributed along trajectories on the surface or inside the Bloch sphere, 
depending on the degeneracy of the spectrum. At certain values of the coin parameters 
the spectrum of the time evolution operator becomes 2-fold degenerate, but there might 
exist unique eigenvalues at the top and bottom of each quasi-energy band. The eigenstates 
corresponding to such eigenvalues are robust against arbitrary temporal variations in the 
bias parameter of the coin, as long as rest of the parameters remain unchanged. 


1 Introduction 

Quantum information and computation is a field which uses concepts from computer 
science to understand the ability of quantum systems to store and transmit information 
and to perform computations. The quantum walk is a prime example of this vast field: 
it was conceived by adopting the concept of random walks from classical information 
theory to quantum dynamics. The quantum mechanical version of the random walks 
shows interference effects which might allow faster algorithms than those based on the 
classical random walks. 

Quantum walks have been studied intensively for over a decade. The discrete-time 
quantum walk was introduced as a unitary process where a single particle, called the 
walker, hops on a lattice and the direction of movement is decided by its internal state, 
called the coin degree of freedom mm- The propagation of the walker is ballistic as op¬ 
posed to classical random walks where the walker propagates diffusively. The proposition 
of quantum walks led to studies on whether this effect could be used to design compu¬ 
tational algorithms with a quantum speedup The simplicity of the quantum 

walk model has also inspired many studies on various aspects of quantum dynamics and 
quantum information, such as localisation mm, quantum to classical transition jljj and 
quantum state transfer HDunmanai. Moreover, it has been shown that the quantum 
walk can be used to perform universal quantum computation dH5]. Overviews on the 
properties of quantum walks can be found in Refs. pL6 . E]. They have also been im¬ 
plemented experimentally with optical lattices [18, 19], trapped ions [[20l |2T] and optical 
photons mmm- 

The quantum walks come in many varieties. Ultimately the walks are defined by 
the underlying lattice structure and the transition rules for moving between the lattice 
sites. Together they define the time evolution operator, which describes the evolution of 
the wave function between subsequent time steps. The behaviour of the walk is entirely 
determined by the eigenvalues and -vectors of the time evolution operator, and the initial 
state of the walker. The purpose of this paper is to explore the variety of the spectrum of 
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the time evolution operator with arbitrary 2-dimensional unitary coins for translationally 
invariant discrete-time quantum walks on the N-cycle. The eigenvalues and eigenvectors 
for this problem have been solved by Tregenna et al [22], who also derived formulas for the 
asymptotic time-averaged probability distributions in terms of arbitrary coin parameters. 
In this paper, we show that by choosing an appropriate parameterization for the coin, 
each parameter plays a distinct role in the properties of the eigenvalues and eigenstates. 
The solution allows therefore to design quantum walks which have precisely the desired 
spectral features. For example, the gaps between the eigenvalues determine important 
properties of the walk, such as mixing time [2J, and it is desirable to have control over 
such features. 

To further investigate the properties of the eigenstates, we obtain the Bloch represen¬ 
tation of the reduced eigenvectors in coin space, and study the properties of the eigen¬ 
vectors as a function of coin parameters. The eigenvectors are found to be distributed 
along well-defined trajectories on the surface or inside the Bloch ball, depending on the 
degeneracy of the spectrum. With special attention to degenerate eigenvalues, we point 
out that there may also exist unique eigenvalues in an otherwise degenerate spectrum. 
The eigenstates corresponding to the unique eigenvalues have very special properties, 
since they remain fixed for arbitrary variations of the bias parameter of the coin. These 
states are thus protected from arbitrary temporal variations of the bias parameter, if the 
rest of the parameters remain fixed. 

This article is organized as follows. The discrete-time quantum walk and the param¬ 
eterization of the coin is introduced in Sec. [2] The eigenvalue spectrum is studied in 
Sec. [3j and the Bloch representation of the eigenvectors is given in Sec. [4] The protected 
eigenstates are discussed in Sec. [5] and the conclusions are presented in Sec. [6] 


2 Discrete-time quantum walk 

The discrete-time quantum walk proceeds by repeated steps of coin flipping and transla¬ 
tion to the left or right, depending on the state of the coin mm- The walk is assumed to 
take place on a 1-dimensional lattice with periodic boundary conditions, and translation 
invariance is imposed by taking the same coin at every site of the lattice. The wave func¬ 
tion of the walker can be written as = 'Y2 xc ipxc(t) I 2 -)! 0 ), wh ere the amplitudes 

-0 xc (t ) are complex numbers. The number of sites is N and the coin variable c can take 
values {0,1}. The total Hilbert space is a tensor product of the spatial and coin spaces, 
|^(t)) € Hx <8>Hc, with dimension 2N. 

The time-independent step operator U evolves the wave function by one time step: 
\ip(t + 1)) = U\ij}(t)'). It can be written as U = TJ-, where the coin flip acts on the 
coin degrees of freedom only: T = / <g) F. The coin flip F is represented by a 2- 
dimensional unitary matrix. Such a matrix has generally 8 real components, since each 
of the 4 complex elements can be represented by 2 real numbers. Requiring the matrix 
to be unitary, FF t = F^F = /, imposes 4 independent conditions for the elements, 
and therefore the matrix is represented by 4 independent real parameters. One of the 
parameters is just the total complex phase of the matrix, and it will be ignored since 
it does not have any effect on the time evolution of the probability distribution of the 
walker. The rest of the parameters can be chosen in infinitely many ways, but as will be 
seen below, the following parameterization is particularly insightful when studying the 
spectrum of U: 

F _j 0 ( VR c ia VT^Re~'P\ 

' V-01 - R VRe~ ia ) 

The range of the parameter R is [0,1] and it is known as the bias parameter, as it sets the 
weight between the coin components 0 and 1. The angular variables are on the interval 
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[0, 2n). The coin flip matrix introduces mixing between the coin components |o) and 1l), 
which correspond to left and right directed motion of the walker, respectively. 

The conditional shift operator T moves the walker to the left or right depending on 
its coin state. It can be written concisely as 

T = T_ ® P 0 + T+ ® Pi (2) 

where the operator P c = |c)(c| is a projector to the coin state c, and T± = 

|x ± l)( ;r | is the translation operator. Periodic boundary conditions are imposed by 
defining the translations at the righthand boundary as — l| = |0)(iV — l| and at 

the lefthand boundary as | — l)(0| = N — l)(0|. 


3 Eigenvalue spectrum 


The eigenvalues of the time evolution operator U can be solved by standard methods of 
quantum mechanics. Translation invariance implies that the wavenumbers k in Fourier 
space are good quantum numbers, which means that the eigenvalue equation for U be¬ 
comes 2x2 block-diagonal, and each block can be solved independently. Since U is 
unitary, its eigenvalues are complex numbers with absolute value 1 and can be written 
as e*\ The phases A can be interpreted as quasi-energies, if the time evolution of the 
wave function is thought to be generated by a time-independent Hamiltonian H. The 
time evolution operator is then given by U = where the time step is chosen as 

At = 1. Indeed, the quasi-energies A are distributed along single-particle energy bands 
that become continuous as the size of the lattice approaches infinity. The eigenvalues 
and eigenvectors for the IV-cycle with a similar parameterization for the coin have been 
previously solved by Tregenna et al |25j , but the solution is summarized here since the 
notation will be useful later. 

The eigenvalue problem for the time evolution operator is given by the equation 

U\^ x ) = e iX \^ x ), (3) 


where the eigenvectors are denoted by |^a) = Ylx *( x > c)|*)|c) = \k)\c) 

and the Fourier-transformed amplitudes are 'ip\(k,c ) = ^2x=o G. l ‘ 27Tkx ^ N 'i/j\(x, c). The 

task is to find the eigenvalue phases A £ [—7T, it) and the components of the eigenvectors 
%j)\. In the Fourier basis, the translation operators and T + act diagonally and the 
wave function is just multiplied by a phase e ~ l2 ' irk / N and e* 27rfc / i?v . respectively. Thus, the 
2 N equations in Eq. ([3j factorize into N pairs of decoupled equations, involving only the 
components ip\(k, 0) and i^\{k, 1) at each wavenumber k: 


f9oo (Mk,0)\ =Q 

\5io 5n/ \^a(M)/ 


where the matrix elements are given by 

5o0 (A, k) = Vr e i(«+0- 2 *£) _ (5) 

g 10 (k) = (6) 

9oi(k) =VT^Re~ i2 ^ (7) 

gu(X,k) - e a (8) 


for fixed coin parameters R, a , /3. Denoting the matrix in Eq. Q as G, this equation 
has non-trivial solutions if and only if the determinant of G is zero. For most choices 
of the coin parameters, non-trivial solutions exist for only one wavenumber k, and the 
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Figure 1: Distribution of the eigenvalue phases A (k,z) as a function k and z. 
a) Deformation of the spectrum as R varies from 0 to 1. The eigenvalues with 
the same R are marked with the same colour; the negative (positive) values 
correspond to z—1 (z= 2). The system size is N=20 and the values of the other 
coin parameters are a—0, /3= 0. b) Dependence on the coin parameters a and (3. 
The system size is iV = 10 and R— 0.5. For the lower bands z—1, /3= 0, and for the 
upper bands z— 2, a=0. 


eigenvectors are strictly localized at each k. There exist however special points in the 
parameter space where the spectrum becomes degenerate and Eq. (jdj allows solutions 
for two different wavenumbers, as discussed below. 

The eigenvalues e* A can be solved from Eq. Q for each k, and there exist two distinct 
solutions labeled by an integer z = {1, 2}: 


g iX(k,z) 


JP 


Vr cos (a 


2nk/N) + i (—l) 2 i/l — Rcos 2 (a — 2 nk/N) 


(9) 


where the term in the square brackets is a complex number with modulus 1. The 2N 
eigenvalues are thus labeled by the wavenumber k = {0,1,..., N — 1} and the integer z. 
The eigenvalue phases A (k,z) are plotted as functions of k and 2 in Fig. [l] for different 
values of the coin parameters. The eigenvalues for different k are distributed on two 
well-defined quasi-energy bands, labeled by the integer z. The shape of the bands is 
entirely determined by the coin variables R, a and (3, and changing the total number of 
sites changes only the relative distances of the eigenvalues along the band. As N —> oo, 
the eigenvalues on each band become infinitesimally close to each other. 

Each of the parameters R, a and (3 correspond to different properties of the spectrum. 
As seen in Fig. 0 *). the parameter R determines the general shape of the spectrum, 
interpolating between a completely flat spectrum at R = 0 and a linearly spaced spectrum 
at R = 1. In the former case the eigenvalues are IV-fold degenerate, while in the latter 
case the gaps between the eigenvalues are degenerate. In the special case R = 0, Eq. © 
implies that the phases of the eigenvalues are given by 

A (k,z) = p + {-l) z | Vfc; R = 0, (10) 

and for R = 1 they become 


A (k,z) = 


/3 + (— l) 2 (a — 2 irk/N) if sin(a — 2nk/N) > 0 
/3 — (— l) 2 (a — 2nk/N) if sin(a — 2t xk/N) < 0 


R = 1. 


( 11 ) 


Eq. (11) shows clearly that the eigenvalue phases for R = 1 are distributed symmetrically 


around the value /3, with equally spaced intervals AA = 2tt/N. 

The role of the parameter (3 is just to set the “base point" of the spectrum, as can be 
directly from Eq. ©• The effect of changing f3 is that each eigenvalue phase is shifted 
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by a constant amount A A = A/3, as illustrated in Fig. 0 >) for the z = 2 eigenvalues. 
The role of the parameter a is slightly more subtle. As shown in Fig. 0 >) for the 
z = 1 eigenvalues, increasing the value of a increases A for one half of the spectrum, and 
decreases A for the other half. Changing the value of a introduces asymmetry between 
the left and right parts of the spectrum. If the change in a is an integer multiple of 
2n/N, this corresponds to a mere relabeling of the wavenumbers k, which is evident from 
Eq. (J9]) . Therefore the spectra corresponding to a = —tt/N and a = n/N in Fig. [ljn) 
are equal, except that the values are shifted horizontally by A k = 1 . 

For most values of a, the eigenvalues on the left and right side are not equal, and the 
spectrum is nondegenerate. Thus, the left-moving and right-moving components of the 
walker are distinguishable, and the walk is called chiral. 


4 Bloch representation for eigenvectors 


In the previous section, it was seen that each coin parameter R, a, (3 has a distinct role in 
the eigenvalue spectrum of U. It turns out that the coin parameters also control different 
features of the eigenstates. We obtain a geometric picture of the eigenstates which proves 
to be very useful in the analysis of the eigenstates. 

The eigenstates take different forms in the degenerate and non-degenerate cases. The 
condition that two eigenvalues are equal, e WM) = e *A(fcV) fo a( fo tq ie following relations: 

cos(a — 2n k/N) = cos(a — 2ixk! /N) (12) 

(-1 ) 2 V 1 - 3?cos 2 (a - 2-Kk/N) = (-1 Y' y/l - Rcos 2 (a - 2nk'/N) (13) 


The second condition is just equal to z = z 1 , since the signs must be equal on both sides. 
The first condition allows two solutions for the wavenumbers, either k' = k or 


k' = —k + a — 
7r 


mod N. 


(14) 


Since the wavenumbers are labeled by integers, the above equation only makes sense if a 
is an integer multiple of tt/N. Therefore the spectrum becomes degenerate at 2 N discrete 
values of a: 

77 n E {0,1,..., 2A- 1}. (15) 


a = n 


N 

A similar result has also been obtained by Tregenna et al [25]. In the maximally degen¬ 
erate case, there are N distinct eigenvalues, which are labeled by N/2 pairs of conjugate 
wavenumbers (k,k') and the integer z, and the eigenvalues are labelled by A (k,k',z). 

Recall from Eq. Q that the components of the eigenvectors for wavenumber k can 
be non-zero iff det(G) = gooSn ~ 510501 = 0. By a straightforward calculation, it can 
be shown that this condition is equal to Eq. (12). Therefore, if the spectrum is non¬ 


degenerate, only the k' = k component of the eigenvectors is non-zero, and the eigenvec¬ 
tors are localized at wavenumber for each eigenvalue A (k,z). If a is an integer multiple 
of tt/N, the AZ-component of the eigenvectors is also non-zero, where the conjugate com¬ 


ponent k' is given by Eq. (14). 


The relations between the two coin components are trivially solved from the first row 
of Eq. Q: 


V’A (k, 1) 
i'xik', 1 ) 


5oo(A, k) 
5oi (k) 


Tp x (k,0) 


5oo(A, A:') ~ , 

7jr~ Tpx(k , 0 ) 


5oi (A' 

c) = 0 V c, k” 7 ^ k, k'. 


(16) 

(17) 

(18) 
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where |< 701 1 > 0 if R < 1 and A = A (k,k',z). The condition © refers to the non¬ 
degenerate case only. The second row of Eq. Q is then always true if det G = 0. 

To gain some insight to the properties of the eigenstates, it is useful to move to a 
geometric picture. The reduced density matrix of the coin, pc = Trx (|V’a)(' ! /’a|) , is a 
2x2 matrix that can be represented by a vector r. pointing to a point inside the Bloch 
ball [26]. The density matrix is given by pc = |(/ + r ■ <r), where the <r,; are the Pauli 
matrices. If the spectrum of U is non-degenerate, the eigenstates are product states 
of wavenumbers and coin states, the reduced coin eigenvectors are pure states, and the 
Bloch vector points always to the surface of the sphere. In the degenerate case, the 
eigenvectors are not necessarily product states, and the Bloch vector might be inside the 
sphere. To obtain an explicit representation for the eigenvectors, it is necessary to obtain 
an ortho normal eigenbasis. 


From Eqs. (16-18), the eigenvectors for degenerate eigenvalues can be constructed as 


IV’a) = s\k) 


Iq\ _ k) | > 

1 ’ 9oi(k) 1 / 


■ e^s'lA:') f|0) - 


flop (A,//) 
goi(k') 


ll> 


(19) 


where s and s' are real numbers and the complex phase of the | fc) |0)-component is set to 
0. The variables s, s' and uj are free parameters, and any choice of these parameters gives 
an eigenvector of U, but s and s' are not independent once the vectors are normalized. 
The normalization (V'aI^/’a) = 1 imposes the constraints 


s = 


\9oi(k' )\ 2 - [|ffoo(fc)| 2 + |goi (fc)| 2 

\goo(k')\ 2 + \goi{k')\' 


( 20 ) 


0 < s < s max 


IffOl^Ql 

V / l5 f oo(A,/c)| 2 + \g 0 i{k)\ 2 


( 21 ) 


where the limit s —> s max corresponds to s' —> 0. Any vector satisfying these relations 
is a normalized eigenvector of U. but two eigenvectors with parameters wi} and 

{s 2 , s' 2 , 0 J 2 ] ar e not necessarily orthogonal to each other. But it is always possible to pick 
two orthogonal eigenvectors for each doubly degenerate eigenvalue, if these parameters 
satisfy 


S2 


UJ2 


\goim 


\goo(k)\ 2 + \g 0 i(k)\ 2 


- st 


= uq + 7T. 


( 22 ) 

(23) 


Any choice of si and uj 1 gives then an orthonormal pair of eigenvectors. Note that the 
weights si can be different for each eigenvalue, as long as they fall on the interval given 
in Eq. (21). 

The Bloch sphere representation of the reduced coin eigenvectors can now be con¬ 
structed from Eq. (19) with the result 


r x 


r z 


= — 2 Re(0) 

= 2 lm(0) 

goo (A, A:) 2 \ 
goi(k) ) 



goo (A, k') 2 \ 
goi(fc') 2 


where the function 0 is dehned as 

q_ 2 goo (A, k) , 2 g 00 (\,k) 

~ s iMkV goi(k') 


(24) 

(25) 

(26) 


(27) 
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Figure 2: Distribution of the reduced coin eigenvectors on the Bloch sphere for 
different values of the coin parameter R. The labeling of the markers is the same 
in all graphs. In graph c), the filled points correspond to z—1 and the unfilled 
points correspond to z—2. The weights between the conjugate components are 
set equal, s\ = S{. so that there is no difference between the z components in a) 
and b). The system size is N —20 and a=0, /3=0. 


The Bloch ball representation of the eigenvectors is plotted for different values of the coin 
parameters in Fig. [2j where we have moved to the spherical coordinates (r x , r y , r z ) —> 
(r,9,(j>). The eigenvectors are distributed along well-defined trajectories, which become 
thicker as the number of lattice points increases, in analogy to the eigenvalue bands 
discussed in the previous section. 

It can be seen from Fig. Ht>) that for small values of R, the coin eigenstates are close 
to the equator where 9 = 7 t/2 , and increasing R brings them closer to the poles of the 
Bloch ball. Fig. [2jc) shows that the eigenstates with the same k but different z are at 
approximately same distance from the point <f> = ir. The z components of the eigenvectors 
traverse around the axis in opposite directions, with the angle cj) increasing or decreasing 
depending on z. This is in contrast to the non-degenerate case, where the angle 4> can be 
shown to be given by <fi = /3 — f + 2-/r^, ie. independent of z. Note that the eigenvalues 
corresponding to k = N/2 ± 5k are equal (see Fig. [I]), and the eigenvectors for the same 
value are distributed at equal distance from the equator. These considerations hold when 
the weights s and s' are set equal in Eq. (19), and it should be noted that changing the 
parameter s can also deform the trajectories, but these deformations are not explored 
here. 

The degeneracy of the spectrum of U is a signal of some kind of redundancy in 
its representation. This redundancy implies that there exists a non-trivial symmetry 
operator S that leaves U unchanged: 


SUSP = u. 


(28) 


7 

















With the solutions for the eigenvectors of U, such a symmetry operator can be constructed 
by considering the orthogonal eigenvectors for each degenerate eigenvalue A (k,k',z). 
Since the operator U is “blind” to the transitions between degenerate eigenstates 
and IV'a)) a non-trivial unitary operator that satisfies Eq. (28) is given by 


5 = 


H (IRKRI + IRKRI) 


(29) 


Ag A. n 


\eA d 


where A n is the set of non-degenerate eigenvalues (if any) and is the set of all de¬ 
generate eigenvalues. When written out explicitly, the operator S consists of transitions 
between the left-moving and right-moving components k and k!, and simultaneous ro¬ 
tations of the coin. It must therefore be understood as a joint symmetry which acts 
simultaneously on the position and coin degrees of freedom. 


5 Protected eigenstates 


Even if the rest of the spectrum is degenerate, in some cases the highest or the lowest 
eigenvalue might be unique, as seen for k = 0 and k = N/2 in Fig. ER Based on the 
discussion in the previous section, the eigenstates for unique and degenerate eigenvalues 
are characteristically different, since the latter are superpositions of two wavenumbers k 
and k'. As seen in Fig. [2j the eigenvectors for the unique eigenvalues correspond to nodes 
of the trajectories in the Bloch ball, and these eigenvectors remain fixed as the value of 
R changes. 

The robustness of the special eigenvectors can be understood using the results from the 


previous section. According to Eqs. (14 15), the wavenumbers for degenerate eigenvalues 


are arranged into pairs (m, n — m) where m = {0,1,..., M}, and (n + l, N — l) where 
l = {1, 2,..., L}. If n is odd, M = [n — l)/2 and every wavenumber in the first group 
has a matching pair. But if n is even, M = n/2 — 1 and the wavenumber k = n/2 does 
not have a pair, and the corresponding eigenvalue is unique. Similarly if IV — n is even, 
L = (N — n)/2 — 1 and the eigenvalue corresponding to k = (N — n)/2 (modulo N) is 
unique. The values of k with a unique eigenvalue are summarized in the following table: 



N even 

N odd 

n even 

k = {n/2, (N — n)/2\ 

k = n/2 

n odd 

- 

CM 

'e' 

I 

II 


For example in Fig. ER n = 0 and N = 20, and the eigenvalues at k = 0 and k = 10 
are unique. 

The widely used Hadamard coin is a curious special case in the parameter space. In 
the parameterization introduced in Eq. ([!]), it corresponds to the values {.R, a, f3} = 
{1/2, 37 t/2 , 7 t/2 }. If N is an integer multiple of 4, this choice implies that the spectrum is 
degenerate and there are 2 unique eigenvalues in the spectrum, since the value n = 3N/2 
is an even integer, see Eq. (15). If IV is an odd integer multiple of 2, the spectrum is 


fully degenerate since n is an odd number. For odd N. the parameter a is not an integer 
multiple of n/N, and the spectrum is always non-degenerate. 

The robustness of the eigenstates is explained by the behaviour of the matrix elements 
of G. Recall from last section that the components of the eigenstates for unique eigenval¬ 


ues must satisfy Eq. (16). If a = rnr/N and k = n/2, then a — 2irk/N = 0 and it is easy 
to show that goo/goi = (— l) 2 e^ Q+ ^ +37r / 2 / If k = (N — n)/2, then a — 2nk/N = it and it 
follows that 500/501 = (—l) z e^ -Q+ ^ +7r//2 / The eigenstates are thus equal superpositions 
of the coin components |0) and |l), with the phase difference determined only by the 
coin parameters a and (3 (and the eigenvalue label z). This property is reflected in Fig. 
|2|u) where the Bloch angle 9 at k = 0, N/2 is equal to 7r/2 for all values of R. 
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The robustness of the special eigenstates might be used to store quantum states and to 
protect them against arbitrary variations of the parameter R. If the eigenvalues z = 1,2 
at wavenumber k are unique and the rest of the spectrum is degenerate, the eigenstates 
|) an d |'0A(fc,2)) remain “frozen” even if the parameter R changes at every time 
step, as long as R is spatially homogeneous and the parameters a and /3 remain fixed. 
Therefore, if the initial state is 

IV’(O)) = x 0 |<y?) + Xi\lp X (k,l)) + 1 Vf\(fc,2) ) (30) 

where | ip} is orthogonal to |V’A(fc,i)) an d the overlaps between the state at time 

t and the eigenstates remain fixed at all times, up to a complex phase: 

(V’A(M)|V’( i )) = e iyi xi (31) 

{■0A(fc,2) |V>(*)> =e iy2 x 2 , (32) 

and the probabilities to measure the states |V’A(fc,i)) and |V ; A(fc, 2 )) are constant. The 
quantum walk can thus act as a memory which is protected from temporal variations of 
R. 

The robustness of the protected states is reminiscent of symmetry-protected topolog¬ 
ical phases, where the edges between different phases may host chiral states inside the 
band gap. Such states are also robust to arbitrary variations of the system parameters, 
as long as the symmetry class of the Hamiltonian does not change. Such states have also 
been found in the so-called split-step quantum walk mm and also in disordered quan¬ 
tum walks [29]. In fact, the topological structure of quantum walks has been found to be 
richer than that of non-interacting many-body systems, and the topological features in 
the quantum walk are not always captured by standard tools used to characterize topo¬ 
logical phases in non-interacting systems m- In the simple case that we have considered 
here, the properties of the eigenstates were found to be very different in the presence 
of symmetry. The topological features in quantum walks are still not completely under¬ 
stood, and it interesting that even the simplest quantum walk can host states which are 
robust against certain variations of the system parameters. 


6 Conclusions 

The spectrum of the time evolution operator of discrete-time quantum walks was analysed 
for a general 2-dimensional unitary coin. The eigenvalues are labelled by two integers k 
and z, where k is the wavenumber in Fourier space and z labels the upper and lower quasi¬ 
energy bands. If the coin is parameterized suitably, each parameter can be associated 
with distinct properties of the spectrum. The bias parameter R controls the dispersion 
of the bands, the angle a changes the symmetry between left-moving and right-moving 
components, and the angle (3 sets just the base level of the quasi-energies. 

For the eigenstates, the geometric picture using the Bloch ball was found to be very 
useful. The reduced eigenvectors in coin space are distributed along well-defined trajec¬ 
tories, which become deformed as the coin parameters are varied. For small values of the 
bias parameter R, the eigenvectors remain close to the equator. Increasing R brings them 
closer to the poles, with all eigenvectors localized at either c = 0 or c = 1 when R = 1. If 
the parameter a is an integer multiple tt/N, the spectrum becomes degenerate, which has 
a significant impact on the eigenstates. In the non-degenerate case, the eigenvectors are 
product states of the Fourier modes and coin states, but in the presence of degeneracy 
they are entangled superpositions of conjugate wavenumbers k and k!. In the latter case 
the reduced coin states are mixed states, and they are not constrained to the surface 
of the Bloch sphere. The representation of the eigenvectors is also not unique in this 
case, and we have identified the parameters si and as a certain kind gauge freedom, 
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which can be chosen freely while the resulting states are still orthonormal eigenstates of 
U. However, the Bloch representation of the eigenvectors is not equal for different si and 

CJl. 

Finally, depending on the evenness of the number of lattices sites N and the difference 
(N — n), where a = nir/N, it has been shown that there may exist unique eigenvalues 
in the degenerate spectrum. These points are located at the top or bottom of the quasi¬ 
energy bands, and the eigenstates corresponding to the unique eigenvalues are robust 
against arbitrary variations in the parameter R. If the initial state of the quantum walk 
has a non-zero overlap with such states, the overlap will remain constant in time, even if 
the parameter R changes between time steps. Since these states remain constant in time, 
they might be useful for storing quantum states, and such a memory would be protected 
from arbitrary temporal variations in the parameter R. 
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